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We study the large A limit of the loop-dependent characteristic functional Z(\) =< exp(iA <f> v ■ 
dx) >, related to the probability density function (PDF) of the circulation around a closed contour c. 
The analysis is carried out in the framework of the Martin-Siggia-Rose field theory formulation of the 
turbulence problem, by means of the saddle-point technique. Axisymmetric instantons, labelled by 
the component a zz of the strain field - a partially annealed variable in our formalism - are obtained 
for a circular loop in the xy plane, with radius defined in the inertial range. Fluctuations of the 
velocity field around the saddle-point solutions are relevant, leading to the lorentzian asymptotic 
behavior Z(X) ~ 1/A 2 . The 0(1/A 4 ) subleading correction and the asymmetry between right and 
left PDF tails due to parity breaking mechanisms are also investigated. 

PACS numbers: 47.27.Gs, 11.15.Kc 



I. INTRODUCTION 



The study of the statistical properties of circulation in fully developed turbulence has been attracting a large deal 
of attention along the last few years [Q-Q. The main motivation relics on the emergent picture of turbulence as a 
phenomenon intrinsically related to the dynamics of vorticity filaments, clearly observed for the first time in direct 
numerical simulations of the Navier-Stokes equations Q. Filamentary structures seem to have, in fact, a fundamental 
place in the hierarchy of eddy fluctuations, as advanced in a recent phenomenological work of She and Leveque [^) , 
where multifractal exponents of velocity structure functions were predicted to very accurate precision. 

An earlier theoretical analysis of the problem of circulation statistics was attempted by Migdal pL who proposed, 
using functional methods originally devised for the investigation of gauge theories, that in the inertial range the 
PDF of the circulation T, P(r), evaluated for a closed contour c, should depend uniquely on the scaling variable 
r/A( 2k ~ 1 ^ 2k , where A is the minimal area enclosed by c and k is an unknown parameter. It was initially thought, 
in order to compute fc, that the central limit theorem could be evoked to regard T as a random gaussian variable 
obtained from the contributions of many independent vortices. Using, then, the definition of circulation, 




and the Kolmogorov scaling law, < \v(x) — v(y)\ >~ \x — y] 1 ^ 3 , a simple guess would be k = 3/2, leading to 

< r n >~ A 2n / 3 . We now know, however, from a numerical analysis by Cao et al. that although there is some 
support to the minimal area conjecture and the general existence of a scaling variable, as defined above, the gaussian 
description of the circulation statistics in the inertial range and the "Kolmogorov" exponent k = 3/2 are both ruled 
out (the numerical results indicate k < 3/2). Gaussianity holds only in the integral scales, while it turns out that for 
loops contained in the inertial range the correlation between vortices cannot be neglected, a fact that obstructs an 
application of the law of large numbers. Intermittency, as found in ref. [2] , is signaled at the tails of the circulation 
PDFs, which are fitted by stretched exponentials like P(r) ~ exp (— /3|r|"), where a ~ 1 in the inertial range and 
a ~ 2 in the integral scales. On the other hand, the circulation PDF cores are gaussian, as one could expect. 

An important conceptual point, raised in the same numerical simulation and related to the determination of k, is 
whether the moments of T are independent or not from the form of velocity correlation functions. In order to find 

< r 2 >, for instance, it may be useless to know the two point correlation function < v a (x, t)vp(x' , t) >, since the 
contour integrations which appear in the definition of the square of the circulation and the average over realizations 
of the random velocity field may not commute. 

Our aim in the present work is to study the problem of circulation statistics in the inertial range through the 
Martin-Siggia-Rose (MSR) technique B. In spite of the many years passed since its advent, only recently interesting 
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results were obtained from the MSR formalism, concerning the computation of intermittency effects in problems like 
turbulence in the Burgers model and in the transport of a passive scalar |t],|8| . The basic tool employed in these works 
is the saddle-point method, where instanton configurations and fluctuations around them are assumed to contribute 
in a significant way to the evaluation of the MSR functional. As we will see, a computation carried along these lines 
will give us non-gaussian tails for the circulation PDF, with stretching exponent a = 1, in reasonable agreement with 
the numerical findings commented above. 

This paper is organized as follows. In section II, the basic elements of our formalism are set. We define the MSR path 
integral expression from which the circulation PDF may be derived, and work out its instanton solutions. The saddle- 
point action is then computed. In section III, we move to the next natural step, which is the study of fluctuations 
around the saddle-point solutions. We find that fluctuations contribute in an essential way to the asymptotic form 
of the MSR functional. In section IV, we investigate subleading corrections to the asymptotic expression, induced 
by small (and gaussian) fluctuations of the circulation. As a result, we establish, for the PDF of the circulation, a 
relation between the width of its gaussian core and the tail decaying parameter (3. In section V, we study the structure 
of asymmetric PDFs, due to parity breaking mechanisms, like turbulence in rotating systems or under the action of 
parity breaking external forces. We comment on our results in section VI, pointing out directions of future research. 
In the appendixes, we discuss in more detail computations which underly some of the results presented in the bulk of 
the paper. 

II. INSTANTONS IN THE MSR APPROACH 

As largely known, inertial range properties of three-dimensional turbulence may be modeled by the stochastic 
Navier-Stokes equations ||] , 

d t v a + vpdpv a = -d a P + vd 2 v a + f a , 

d a v a = , (2.1) 
where the a — 1, 2, 3 and the gaussian random force f a {x, t) is defined by 

< f a (x,t) >= , 

< f a {S,t)f f3 {x , ,t') >= D aP (x-x')5(t - = £> exp (- |x ~J 1 ) S a0 8(t - t') . (2.2) 

Above, L is the typical correlation length of the energy pumping mechanisms, acting at large scales. The other 
important length in the problem, according to Kolmogorov theory JlQ] , is r\ ~ t^ 3 ^ 4 — > 0, the microscopic scale where 
viscosity effects come into play. 

From the stochastic Navier-Stokes equations one may try to obtain, in principle, any velocity correlation function. 
We are particularly interested to study the characteristic functional 

Z(X) =< cxp(iAr) > , (2.3) 



where T is the circulation evaluated at time t=0, as given by (1.1). The contour c used in the definition of T is taken 
here to be the circumference x 2 + y 2 = R 2 , with z = 0, oriented in the counterclockwise direction. A basic condition 
in our analysis is that R is a length contained in the inertial range, that is, r\ R L. The PDF for the circulation 
may be written from the loop functional as 

1 f°° 

p(r) = — / d\cxp(-ixr)z(x) . (2.4) 

It is appropriate, for the computations which will follow, to consider the analytical mapping A — > — iX in the RHS of 



(2.3). At a later stage we will get back to the original definition of A. The MSR formalism M allows us to write the 



path-integral expression 

Z(X) = J DvDvDPDQexp(-S) , (2.5) 



where 
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S = — i J d 3 xdt [v a (dtV a + vpdpv a — vd 2 v a + d a P) + Qd a Vaj 

+ -J dtd 3 xd 3 x'v a (x,t)D al3 (x - x")vp(x" ,t) - XT . (2.6) 

The MSR technique may be used to derive, in an alternative way, the Wyld diagrammatic expansion Ej for 
the computation of correlation functions of the velocity field, obtained directly from the stochastic equations (2.l|). 
This expansion is constructed by taking the non-linear term in the Navier-Stokes equations, related to convection, 
as a perturbation. For this reason, the perturbative MSR- Wyld approach has been criticized along the years, as an 
inappropriate tool to deal with the singular configurations of the velocity field, which are fundamental in turbulence. 
However, the advantage of the MSR formalism is that non-perturbative issues may be addressed in principle, if one 
knows how to find specific configurations of the flow that could represent relevant contributions to the functional 
integration for Z(X). This is precisely the task to which the saddle-point method is devised for. 

The role of the P and Q fields in the above path-integral summation is just to assure that d a v a = d a v a = 0. These 
incompressibility conditions are in fact two of the four saddle-point equations obtained from the action S, viz., 

H = d a v a = , (2.7) 

H = d a v a = . (2.8) 

The other two saddle-point equations are given by 

s s (5r 

j— = i (d t v a - VfjdaV/3 + vpdpia + vd 2 v a + d a Q) - Xj— = , (2.9) 



SS 



i (d t v a + vpdpva - ud 2 v a + d a P) - / d 3 x"D af3 (\x - x'\)vp(x' , t) = . (2.10) 



We have 



<5r s f 

— = — f vpix", Q)dx' p = e 3 /3a — S(r± - R)6(z)S(t) , (2.11) 

0V a OV a J c T± 

where r± — (x 2 +y 2 ) 1 ^ 2 ■ The importance of the saddle-point equations is that they provide a systematic way to study 
the large A limit of Z{X). However, the saddle-point action computed in this way necessarily depends on A in a way 
incompatible with observational results Q). In order to understand it, we observe that the saddle-point equations are 
invariant under the scaling transformations hv a , P hP, Q -> h 3 l 2 Q 

and A — => hX. These relations imply that the saddle-point action has the general form = A 3 / 2 f{X^ 1 ^ 2 v). Since 
we expect to have finite answers in the limit of vanishing viscosity, it follows that ~ X 3 / 2 . This dependence on 
A is exactly the one found in Burgers turbulence for the statistics of velocity differences 0^2| , which we know not 
to reproduce, even qualitatively, the PDFs of the circulation in three-dimensions. A similar difficulty was in fact 
noticed in the investigation of velocity structure functions in incompressible turbulence by means of the saddle-point 
method ||. In order to find physically meaningful PDF tails of the circulation, a solution of this problem will be 
pursued here, based on the definition of an additional field in the MSR path-integral, parametrizing an infinite family 
of saddle-point configurations. 

We would be tempte d to study the above saddle-point equa tion s by fi rst el iminating the P and Q fields in (|1]) and 
( [2.10 ) with the help of (2.7) and ( |2.8| ). All nonlinear terms in ( |2.9| ) and (2.10) would consequently appear projected on 



transverse modes through the use of the tensor Ii a p = d~ 2 (d a d/3 — 8 a p)- However, this is not an adequate procedure to 
follow, in view of the simplifications inherent in the implementation of the saddle-point method to the MSR formalism. 
The central point is that we will be dealing with linear approximations for the velocity field, as a consequence of the 
small radius R of the contour c, in comparison with the large scale length L. We have, thus, 

v a (x,t) =<r a p(t)xp , (2-12) 

with a aa = (due to d a v a — 0). Coordinate independent terms are not written above, since we may impose, from 
invariance under the grou p of time-dependent translations, the saddle-point solution to satisfy v a (x = 0, t) = Q (see 



appendix A). Using ( 2.1 2| ) we observe that expressions like Hapv-fd^vp, related to the global nature of the flow, would 



not be precisely defined. A simple way out of this problem, usual in applied mathematical studies of the Navier-Stokes 
equations jll| , is to write the pressure as a quadratic form, 
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P = A aj3 x a xp , (2.13) 

so that d a P exactly cancels in (2.10) any s ymmetric tensor acting on the spatial coordinates, which would appear in 
the linear approximation. Therefore, ( 2.10|) may be written as an equation for the time evolution of the antisymmetric 



part of the strain field, 

d f 

—a s a/3 + (a s <x s + a s a s ) al3 - i j d 3 xd [f3 D a]l (\x \)v y (x, t) = , (2.14) 

where we have defined 

a a/3 = -^(VasP + Vf3a) , (2-15) 

a i/3 = l^PafS ~ a P<x) > ( 2 - 16 ) 

d [0 D ah (\x\) = 1 (dpD ai (\x\) - d a D^QS\)) . (2.17) 

An important remark is that ( ^.12| ) is not assumed to represent a direct modeling of the velocity field in sustained 
turbulence, which we know to be associated with many different length scales and singular structures. The idea of 
the instanton method, as advanced by Falkovich et al. J8|, is in fact to consider, in the MSR framework, smooth 
configurations and perturbations around them that may condense some information on the statistics of the strong 
(intermittent) fluctuations of the velocity field. The situation here is analogous to the well-known instanton approach 
to the double well potential in quantum mechanics ||l4j| , where instantons are obtained as saddle-point solutions, 
yielding extremes of the euclidean action. It is clear in that case that the smooth kink/anti-kink form of the instanton 
configurations cannot be taken as a direct representation of the quantum-mechanical dynamics, which has a picture 
as a sum over particle paths with complex weights exp(iS'). In the turbulence context, instead of transforming time 
into an imaginary variable as it is done quantum mechanics, we look for saddle-point solutions, considering, in the 
MSR action, the analytical mapping A — > — iX. A deeper analogy, which should also be noted, is provided by the 
phenomenon of localization in condensed matter physics. There is, in this case, a functional integral formalism, where 
smooth instantons may be found, giving expressions for the tails of the density of electron states Jl5[ . The similarity 
with the turbulence problem is a strong one: while in the condensed matter system localized wavefunctions define 
some multifractal set, the same phenomenon takes place in turbulence, regarding the fluctuations of the velocity field. 
Also, the limitations of the instanton method are exactly the same in both problems. Either in localization or in 
turbulence the core of the density of states or of the PDFs, respectively, cannot be obtained from the saddle-point 
technique. To understand it in our analysis of the statistics of circulation, we note that for large values of A the 
functional Z(\) gets its more relevant contributions from the tails of the circulation PDF. At the core, where the 
PDF is essentially stationary, fluctuations of exp(iAr) will tend to prod uce de stru ctive i nterfe rence 



Our problem has been reduced so far to a n analysis of equations (|2.8|) , ( |2.9| ) and ( 2.14 ), where in the second 



equation the velocity field is given by (2.12). Since these equations are invariant under rotations around the z 
axis, it is interesting to look for axisymmetric solutions. In the linear approximation, the most general form of an 
axisymmetric strain field is given by 



a(t) = 



a(t) b(t) 
-b(t) a(t) 
-2a(t) 



(2.18) 



The above form of a(t) has a simple hydrodynamical interpretation. Taking a > 0, for instance, streamlines are 
just expanding spirals which approach in an exponential way the xy plane from both regions z > and z < 0. It 
is important to note that a zz {i) = —2a (t), which has the dimensions of the inverse of time, plays the role of an 



arbitrary external function in eq. (2.14). In other words, vorticity is controlled by stretch, associated to a(t). We 
should try to find instantons (the solutions of the saddle-point equations) for any well-behaved function a(t) (with 
a(t) — > as \t\ — > oo) and then sum up their contributions in the path-integral expression for Z(X). This suggests an 
alternative strategy of computation, where a(t), or some variable related to it, would appear from the very start in 
the MSR formalism as a field labelling families of velocity configurations. There are, in fact, many different way s to 



implement this idea, distinguished essentially by computational convenience. Our choice consists in writing (2.5), up 
to a normalization factor, as 

Z{\) = J DvDvDPDQDa s 6 [d a v p \ z=0 + dpv a \ z=0 - 2a s a/3 ] exp(-S) 

Da s \ DvDvDPDQDQexp(-S) , (2.19) 
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where a s a n = a s a Jx,y,t) and Q a p = Q a p(x,y,t) are symmetric matrices and 



S = S 



- J dxdydtQ aP (x,y,t) ■ (d a vp\ z=Q + dpV a \ z=Q - 2a s a/3 (x,y,t)) 



(2.20) 



The meaning of (2. IE) is that we sum up the contributions to the path-integral expression in two steps: first by 



considering velocity configurations which satisfy d a vp 



z=0 



dnvJ 



2a s a p{x,y,t), for a given field a s aj3 . The 



summation over the fields a s a o is performed afterwards. The linear approximation for the velocity field corresponds, 
thus, to fields a s a g with slow dependence on the x and y coordinates, within the length scale of the order of R, 
while axial symmetry, a condition related to large values of A, is imposed here as a restriction on the configurations 
for cr^g(t). More precisely, we will consider the sum in ( 2.f 9| ) as carried over the space of axisymmetric fields 
a ap(t) = (<W ~ 3<5 a 3<5/33)a(£), in accordance with ( 2.18 ). This corresponds to replacing J Da s (t) by / Da(t) in ( 2.19 ). 
However, this constraint has to be applied with care, since its meaning is linked to configurations of the velocity field 
defined at length scales larger than the loop's radius R. To state it in a different way, the velocity field that enters in 
the above delta functional is in fact a "smeared" field, given by the contributions of wavenumbers k < 

The saddle-point method is to be used in the first step of computation (where a s a(j is fixed) involving the action S 
rather than S. The only modification of the previous saddle-point equations (2.7 - 2.10), as may be readily seen from 
5, is on (2.9), which must be replaced now by 



SS 

Sv a 



(d t v a - vpd a vp + vpdptia + vd 2 v a + d a Q + dp(5(z)Qp a fj ~ A^j- = 



Sv a 



We also have an additional equation, associated to variations of the field Q a p, 

5S 



6Q 



a/3 



= ~i {d a vp\ z=0 + d p v a \ z=0 - 2a s a/3 (t)) = 



(2.21) 



(2.22) 



This equation, however, is before hand solved by ( 2.12 ) and ( 2.18 ). Using ( 2.11 ), ( 2.12 ) and taking the limit of 
vanishing viscosity, we may write (2.21) as 



d t v a - crpaip + ap-fXydpia + d a Q + dp{6{z)Qp a ) 
R)S(z)S{t) . 



«Ae 3ct/ 3 — d(n 
r± 



(2.23) 



We have, therefore, a closed system of coupled equations, given by (U),([0J) and ( HI ' :t 

is important to state 

t he b oundary conditions that the solutions of these equations have to satisfy. Since the viscosity term appears in 
(2.21) with the opposite sign, compared to the one in the Navier-Stokes eq uatio ns, we impose, in order to avoid an 
unbounded growing of the field v a (x, t), that v a (%, t > 0) = 0. In this way, ( f|f ) leads us to 



iXe 



30a El 5(r ± - R)S(z) 



Also, we require that v a (x,t) — > as t — > — oo. The equation for v a (x,t) may be solved through the ansatz 

oo 

v a (x,t) — e 3 f3 a xp5(z) en(t)r"~V")(rj_ - R) , 



(2.24) 



(2.25) 



Tl=0 



where S^(r± — R) = d n 6(r± — R)/dr\. The boundary condition (2.24) reads now 

co(0~) = £A , 

c„(0") = 0, for n > . 



(2.26) 



We find, substituting (J2.25D in (2.22) 



—c + ac Q = , 
at 

— c n + a(n + l)c n + ac„_i = 0, for n > 



(2.27) 
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and Q a p = (8 a p - S a s6ps)Q, with (below, a = 1, 2) 

oo 

d a Q = -2b(t)x a c n (t)r"- 1 S ( - n \r ± ~ R) , 

n=0 

= o . 

The infinite set of equations ( [2.27j ) as well as (2.28) are solved respectively by 
Cn[t) = — e Jo M e Jo 1 J > 

Q(r±,t) = -2b(t) y)c(t) / ^r<5 (n) (^--R) = -2iA6(t)0(r ± - .ReJo dt ' a(t,) ) , 
n=0 ^ 



(2.28) 
(2.29) 

(2.30) 
(2.31) 



where = (1 + \x\/x)/2 is the step function. Taking ( 2.30 ), the infinite summation in ( 2.25 ) may be exactly 
performed. We find the compact result for t < 0, 



v a (x,t) = i\e 3l 3 a — 8{r±_ - Re$o dt a( * ' > )S(z) 



(2.32) 



In order to get some intuition on the singularity in the above expression, we just recall that the quadratic term for 
v a (x, t) in the MSR action is obtained from 

< exp(i J d 3 xdtv a (x,t)f a (x,t j) >/ , (2.33) 

where the brackets denote an average over realizations of the stochastic force field f a (x, t) . Substituting in this average 
v a (x,t) by the saddle-point solution (2.32), we find 



d 3 xdtv a (x,t)f a (x,t) ~ J dt <j> dx a f a (x,t) 



(2.34) 



where the loop integral is taken around the circumference of radius r±_ = Rcxp(J Q dt'a{t')). We see that ( |]3g ) is m 
fact non-vanishing for configurations of the force field that may produce some circulation around the loop r±_ = R, at 
t = 0, through convective proces ses in the fluid. 

Let us consider now equation ( [2.14 ) for the velocity field, which, using the strain field ( 2.18| ), may be written as 

(2.35) 



(2.36) 



b + 2ab + i J d 3 xd [1 D 2 ] 1 {\x\)v 1 {x,t) =0 . 
Substituting the solution for v 7 (x, t) , in the above expression, we obtain 



b + 2ab = -2ttD q X 



'■Io dt ' a ^d(-t) 



In order to have well-behaved solutions for t — * —00, we see, from (2.36), that it is necessary to have in this limit 
L dt'a(t') — > —00. Motivated by the general idea of a gradient expansion, we will restrict our study, as a first 
approximation, to the effect s of time independent configurations given by a(t) = a > 0. Correspondingly, in the 
definition of Z(X), eq. (2. IE), we will have 



Da s 



da 



(2.37) 



A possible physical interpretation of the above replacement is related to the experimental observation of circulation 
as a more interm ittent random variable than longitudinal velocity differences Thus, in the decomposition of 
the strain tensor ( [2.18 ) into symmetric and antisymmetric parts, the latter is actually the quantity which fluctuates 
more strongly in the "background" defined by the partially annealed field a(t). It is worth observing this kind of 
interpretation is usual in a large variety of systems characterized by different time scales, like spin glasses, for instance, 
in the situation where the dynamics of spin couplings is slow - but not negligible - when compared to the typical 
time for spins to reach thermal equilibrium . 
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Equation ( [2.36 ) may be easily solved, yielding 



b(t) = 



nDnX (R 



L 



2a 



-2a\t\ 



(2.38) 



As it could be anticipated, we see that (2.38) represents the well-known phenomenon of vorticity amplification by 
vortex stretching, controlled by the parameter a. Although viscosity does not enter in this expression, vortex stretching 
is bounded, which would not occur in an inviscid flow. The explanation for this behavior of th e ins tanton solution 
follows from the fact that viscosity has been taken into account in an implicit way, through eq. ( [2.24| ), which defines 
V a {0, t) at the initial time t = 0, so that the saddle-point solutions vanish as t — > ±oo. The peculiar property of ( 2. 38] ) 
that will be important in our subsequent considerations is just the factor A/a, relating A and the vortex stretching 
parameter a to the amplitude of b(t). 

The saddle-point solutions we have found for v a (x, t) and v a (x, t) may be substituted now in the action S to give 



n 2 D R 4 A 2 
a 



2L 2 



We note that a straight application of this result would lead to 

7r 2 DoR A A 2 



Z(X) 



da exp( 



2L 2 



(2.39) 



(2.40) 



which is divergent as the integration region extends to a — » oo (above, A has been substituted by iX). This "ultravi- 
olet" divergence is in fact naturally regularized when we also take into account fluctuations around the saddle-point 
solutions, as shown next. 



III. ANALYSIS OF FLUCTUATIONS 



Denoting the saddle-point fields and fluctuations around them by the indexes "(0)" and "(1)", respectively, we 
write 

v a (x, t) = {S, t) + (f , t) , v a (x, t) = (x, t) + (x, t) , 
P(x , t) = p(°) (x , t) + P« (x, t) , Q(x , t) = Q<°) (2, t) + (x, t) , 



Q a p{x,y,t) = Q { ^(x,y,t) + Q K ^(x,y,t) 



(3.1) 



The action is expressed as S = + where S^ ' is given by (2.39), and we have, up to second order in the 
perturbations, 



S« = - i I d 3 xdt [v£\d t vW + vfdpvM + vfdpvW - vd 2 v^ + d a PW) 
+ i(o)( v ^d fj vW)+Q^d a ««] -U dxdydtQ^ (d a vf + d p v^ 



z=0 



+ i / dtd 3 xd 3 x'v^\x,t)D af3 (x-x')v ( l3 1) (x',t) 



(3.2) 



We included in ( |3.2| ), for the sake of completeness, the viscosity term, which in fact will be assumed to vanish in the 
next computations (nevertheless, we have to keep in mind that viscosity, as discussed before, plays an important role 
in the choice of the boundary condition for ia (x, t) at t = 0). 
The integrations over Q^> and Q^l hnply that 



a) d a vi 1) (x,t)=0 , 

b) d a vg\x,t) = 0, 

c) ^d a vf\S,t) + d v^(x,t) 



. 



z=0 



(3.3) 
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If perturbations are written in a form which satisfy these relations, as we will do, then the fields and Q^l 

may be taken out from S^ 1 ). We are interested to find expressions for {x, t) and t)« (x, t), which describe effective 
degrees of freedom. 



The singularity of v„ (x,t) at r± — Re at and z = 0, given by (2.32), represents a ring that shrinks to a point as 

t — » — oo. One could imagine local fluctuations around ia (x, i) given by variations of the vector field defined on the 
ring, 



(x,t) = Va(0,t)5(r-L - i?e at )<5(z) 



where is the azimuthal angle in cylindrical coordinates. Since ip a (6,t) 
series (p a (0, t) = J2n=-oo 



fo 



(3.4) 

+ 27r, t), we may write the Fourier 



, t) = X)^L-oo <P°^ (*) ex p(* n ^)- The incompressibility condition (^b), however, implies that (p^ = 0, for 



and ifa \t) — 5c(t)e^ a fjXp/rj_. Therefore, we are only allowed to consider amplitude fluctuations as 

X/3 



(x,i) 



5c(t)e 3a(3 ^5(r x - Re at )5{z) 
r± 



(3.5) 



An important remark is that the above expression is valid exclusively for negative times, since Va(x,t > 0) = 0. 

We could also take into account perturbations of the ring that would deform its shape, but a little reflection shows 
they may be neglected. Consider, for instance, perturbations of the ring in the xy plane, given by a field T](8,t): 



v a (x, t) = <p a (r±, 6; n)8(r± - Re at + 7/(0, t))6(z) 



(3.6) 



where the above amplitude tp a is a functional of 77 ( 
order in r](9,t) we may write 



1,t) and satisfies to ip a (r±,9;r] = 0) = iXezafjXp/r^. Up to first 



(x, t) 



d9'r)(6',t) 



¥> a (rx,0;v = 0) 



S(r± - Re at )S(z) + tXe^^O, t)S^ (r x 



Re at )6(z) . (3.7) 



The first term in the RHS of thi s eq uation may be absorbed by fluctuations given by (p~i|). Regarding the second 
term, the same steps that led to ( |3.5| ) give us now dgrj(6, t) = 0, that is, the ring is deformed in the xy plane through 
uniform radius variations. It is clear, due to the derivative of the delta function in (3.7) that (3.5) is in fact a more 
relevant contribution at lower wavenumbers. The same reasoning may be extended to generic perturbations of the 
ring's shape. The approximation of neglecting deformations of the ring would be inconsistent if there were small scale 
fluctuations of the veloci ty f ield taking place in a neighborhood of the ring, as we would conclude from the coupling 
of type vv in the action (3.2). However, as it will be shown in a moment, small scale fluctuations of the velocity field 



are contained only in some small compact region surrounding the origin. 

In view of the action of random forces at large length scales (k < L -1 , in Fourier space), we keep, as a first 
approximation, the linear dependence of the velocity field on the spatial coordinates, introducing fluctuations of the 
strain field as 



(x, t) = a a (t) + u}p(t)e a ^x n 



(3.8) 



This linear expression is the only one compatible with the constraints (3.3a) and (|3.3| c) 



If we take Sc(t) — constant, it is not difficult to see, substituting (3__5) and (3.6) in (3.2), that will not depend 



on a a (t) or for t < 0. In other words, we have defined a "zero mode" configuration, which would render the MSR 
path-integral completely independent on large scale fluctuations of the velocity field. The solution of this problem 
consists in considering generic time-dependent variations 8c(t), precisely as we are doing, in accordance with the usual 
procedure fo r the treat men t of zero modes associated to instantons (l4| . 

Relations ( |3.5[ ) and ( |3.8| ) were both defined through arguments based on the assumption that fluctuations around 
the saddle-point have to be local. We observe, however, that they do not exhaust, in principle, the effective form 
of perturbations, which may occur also at smaller length scales. In order to achieve full expressions for ia\x,t) 
and Va(x*,t), it is necessary to take a closer look at fluctuations associated to the dynamics of the action . 
Disregarding the coupling Va\v^ dpVa ) - a self-consistent approximation, as we will see - one may note that fiK 1 ), 
which governs the random behavior of Va(x,t), is the MSR field theory obtained from the stochastic equations 

dtv^ + 4° V« + i# Vg» = vd\^ - d a pM + (3.9) 
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and the constraints (3.3a) and (|3.3|c). The random force fg^ (x, t) , lik e f a (x,t), is defined by (2.2). A criterion to find 



the region of space where small scale fluctuations determined by ( p.9[ ) may effectively occur is based on an analysis of 
the local power supplied to the fluid by the pressure and external forces. In the absence of perturbations, the laminar 
flow is described by the velocity field Va (x, t) , with power density 



V = vW(x,t) (-d a P^(x,t)+ix p J d 3 x'd [f3 p ah (\x"\)v^(x",t) 



(a 2 + 3b(t) 2 )ar 2 ± - 8a V 



(3.10) 



where b(t) is given by ( [2.38| ) and P^ is obtained according to the discussion which leads to eq. ( 2.14 ). Taking 
a > (DqA) 1 / 2 , the b(t) 2 term in the above equation may be neglected. We get 



3a 3 z 2 



(3.11) 



The lower bound (D^ X) 1 ^ 2 for a does not modify the asymptotic form of Z(X). We may check it by considering any 
regularized version of ( 2.4C ), assuming its measure of integration is still dominated by the exponential factor as a — > 0. 
A more physical view on the lower bound for a, which will become clear later, is that in order to evaluate the MSR 
functional Z(X), it is enough to take into ac count saddle-point configurations which have support in the time interval 
At < (D X)~ 1/2 , so that the power density ( |3 . 1 0[ ) turns out to be dominated by the symmetric part of the strain field. 
The extra supply of power density provided by the pressure P^ and the stochastic force fcP is 



Vi=<vP(-a a i*v + fV) > . 

Since the equations and constraints for Va (a?, t) are linear, they are invariant under the substitutions 
v£\x,t) - Dl /2 vW(x,t) , pW(x,t) - D^P^^t) , fW(x,t) - D^fW&t) . 



(3.12) 



(3.13) 



The factor D$ which app ears in the two-point correlation function of the random force fa is now replaced by unity. 
From eq. (3.9) and (3.12) we get, taking v — > 0, 



At x — we obtain 



Da l Vi{$ - 0) - (S aP - 35 a3 Sp 3 )a < v^(0)vf\o) > 



> 



(3.14) 



(3.15) 



Since a > (DqX) 1 ^ 2 , we have Va\x,t) ~ a(x a — 3S a3 z), which means that the stochastic equation (3.9) involves 
essentially only two dimensional parameters: a and L. Through simple dimensional analysis we may write 



< 



vV{x)vf{x)>=^, 



(3.16) 



where C Q( g is a dimensionless constant. We find, from (3.15) and (3.16), 

Pi = cD , 



(3.17) 



where c = C\\ + C22 — 2C33. From rotation symmetry around the z axis, we have C\\ = C22, and consequently 
c = 2(Cn — C33). Due to the strong anisotropy in the system described by (3.9), we expect to have c ^ 0. 

Con sidering now |x| 7^ 0, we may use dimensional analysis once more to write for the first term in the RHS of 

(Hi, 



\vfdp{< v^vW >) „ {Xa - ZS a3 z) T 



(3.18) 



where C a is a dime nsion less constant. Thus, for \x\ <C L, the RHS of (3.14) is still dominated by the secon d ter m, 
leading us again to ( |3.17 ). It is important to observe that in the analysis presented above, the derivative in ( 3. 18] ) is 
assumed to be a smooth function of the spatial coordinates, a condition that may be not valid in some specific set of 
points, as in a vortex sheet. 







Eq. ( |3. 17 ) is in fact a result similar to the one that would be obtained from a loose application of Novikov's theorem 
7j. We expect stronger fluctuations of the velocity field for positions where \Vo\ < that is fis|| , 



\a 6 rj_ -8aV| < \c\D . 



(3.19) 



The above inequality is satisfied in a region of space bounded by three disjoint surfaces generated by the revolution of 
hiperbolae, as shown in fig. 1 . It is consistent to assume the surfaces have a well-defined meaning only at length scales 
contained in the inertial range. Since a > (DoX) 1 / 2 , we can see that for large enough values of A, the surfaces enclose 
some region fi surrounding the origin, with typical size Ro ~ (\c\Dq/ a 3 ) 1 / 2 <C R- The condition on A is given by 



A 3 A) 



«i? 4 



(3.20) 



This relation defines, therefore, what is meant by the "large A asymptotic limit". 

To construct an effective picture out of these considerations, we imagine that in additional fluctuations of u£ (x, t) 
and Va\x, t) are superimposed to the previous expressions (3.5) and (|3.8|). Physical results are then obtained in the 



Ro/R — > limit. In practical terms, this amounts to rewriting in a form which explicitly takes into account the 
length scales involved here, Rq and R. With this aim in mind, it is useful to employ the following notation: 



v<(x,t), iff^O 
(x, t) , otherwise. 



(3.21) 



Analogous definitions are provided for v2 (x,t). We get, from ( 3.21 ) and (|3 



d 3 xdtv> (d t v> + vfdftv> + vppvW) - i I d 3 xdtv< (d t v, 



Sen 



(0)) 



S01 



+ vfd () v< +v<d p vW) + \ ( dtd 3 xd 3 x«v>(x,t)D a p(x-x«)v>(x',t) 



u 
1 

2 ,/~ ™/ f 



2 Jx,x'GQ 



dtd 3 xd 3 x'v< (x, t)D a p(x - x')v^(x', t) + 



dtd 3 xd 6 x'v>{x,t) 



Sen, 3' 02 



x D afj (x - x')v<(x',t) 



(3.22) 



Acco rding to the above discussion, we take now v^(x,t) and v^(x,t) to be given b y th e former expressions (3.5) and 
(3.8), respectively. On the other hand, at smaller length scales, given by |x| < Rq, (3^) is not expected to reproduce 

the behavior of (x,t) anymore, so that another parametrization is needed, viz., 

v>(x,t) = d a (t) +b a p(t)x[3 . 



(3.23) 



The linear expressions for Q and (|3~2^ ) are associated to the fact tha t we are considering velocity fluctuati ons to 
depend essentially on wavenumbers given by k < L^ 1 and k ~ Rq 1 . Eq. ( |3.23 ) is not constrained by condition (3.3c), 
since it describes fluctuations at length scales Ro <C R. The surface 9f2 which encloses may be viewed as a vortex 
sheet for the velocity field Va\x,t). In appendix B, it is shown that Q is necessarily a sphere of radius Ro, whereas 
b a /3{t) is an antisymmetric tensor and a a {t) — d a (t ) . As the coordinate independent field a a {t) (= d a {t)) may be 
absorbed by pressure flucutations in the action (|3.2|), we may take 



v<(£,t) = u>p{t)e a07 x^ 
v a(x>t) = <l>0(t)e a 7 x. 



7 ' 



(3.24) 



where it)p(t) and <fip(t) are proportional to the vorticity outside and inside f2, respectively. At this point we note that 
(|2.32[) and (|3 24j) give 



d 3 xv£HvWdpvW) = f d 3 xv^(v<d p v<) = 



(3.25) 



proving the self-consistency of the simplification discussed before eq. (3.9). 

From (3.24) we see that (x, t) gives no stretch. This peculiar re sult is relat ed t o the fact that velocity fluctuations 



at length scales larger than R have to satisfy both the constraints (|3.3|a) and (3.2c), which makes the flow described 
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by (3.9) somewhat unusual, when compared to the ones commonly mode led i n fluid dynamics, where condition ( |3.3| c) 
is not imposed. On a more physical ground, we may say the constraint (|3.3| c) means that the symmetric part of the 
strain field is "frozen" and does not fluctuate around the saddle-point solution, which is a natural assump tion, since 
we ta ke it to represent the slow degrees of freedom. We also note that there is no contradiction between ( 3.17 ) and 
(3.24), since a coordinate independent field, as commented before, is not written explicitly for v<(x,t) and v^(x,t). 

We found expressions for v^(x,t), v^(x,t) and v^(x,t), but nothing was said about v^(x,t). As a matter of fact, 
this field will be rep laced, as shown be low, by li near combination of its moments c a p{i) = J dxv^{x,t)xp. 

Substituting Q (= v<{x,t)) and <&2% in d|||), we find, after a lengthy and straightforward computation, 



Z(X) 



(£>oA)i/2 



da / D[Sc{t)]D\p(t)] D[c a {t)]D[4> a {t)\ 



n 2 D Q R* A 2 



dt[c 3 (t)(Mt) + 2afa(t)) + ci(i)(0i(t) - o0i (t)) 

r°° 4 

+ca(t)(0a(i) - afa(t)) + iiR 2 Sc(t)(p(t) + 2ap{t))] - D \ dt[—{c\(t) + c 2 (t) 



^ 2 {t)) + 2, 2 R 2 {^f8c 2 {t)]} 



(3.26) 



where 



c a (t) = ^£a/37 / d xi>p (x,t)x~f , for a = 1, 2 , 



c 3 (t) = nR 2 Sc(t) + \ [ d 3 x[v{-( 
2 J sen 



x,t)x2 — v^(x, t)x\ 



p{t) = uj 3 (t) - <f> 3 (t) . (3.27) 

A simplifying prescription has been used to get fl3.26j ). The exponential factor exp(at) has been removed from the 
expression for (x, t) and the time integrals have been defined for — oo < t < oo. The point in doing so is that we get 
gaussian integrals over Sc(t) and c Q (i), which may be exactly computed. The only consequence of this approximation 
is just a slight and unimportant deviation for the values of coupling constants. Taking into account the boundary 
conditions p(±oo) — Q (±oo) = in the resulting path-integral, the time variable is th en rest ricte d to —1/a <t<0, 
where the saddle-point method is assumed to work (this follows naturally from (2.32) and (2.38), which show that 



-,(0) 



(x, t) and b(t) have lifetimes of the order of 1/a and l/(2a), respectively). We will have, therefore, 



Z(X) 



p OO f- 3 

' da / D[p(t)] TTD[0 a (t)]exp{- 

(DoA)i/2 J « =1 



n 2 D R 
2L 2 

L , 



a 



L 2 

Jd~ 



1/a 



dt[<j>j(t) 



+ 4a 2 ^(t) + <t>((t) + a 2 4>i(t) + <f> 2 2 (t) + a 2 ^) + 2{-)\p 2 (t) + Aa 2 p 2 {€))}} , 



(3.28) 



an expression which involves a set of uncoupled on e-dim ensional harmonic oscillators with coordinates <fi±, <f>2, (j>3 an d 
p. Observe that oj\(t) and u>2(t) do not appear in (3.28 ) . This means that at length scales of the order of R, velocity 
fluctuations are essentially axisymmetric. As smaller length scales (of the order o f Rp ) are considered in the action, 
vorticity fluctuations in all directions of space become important. We may write (3.28) as 



Z(X) 



x G({cf)2\(i)2};-,a 
a 



f oo f 3 

' da I dpdp I I d(f) a d(f) a exp( 

(Da A) 1/3 J £Ji 



tt 2 D q R 4 
2LP~ 



-)G({0 1 |0 1 };i, 
a a 



2D, 



^-)G({0 3 |0 3 };-,2a, 
ZDq a 



-^)G({p|p};i,2a,^(^) 4 ) 
ZL)q a Dq R 



where 



G{{x 2 \x 1 };T,uj,m) 



2?r sinh(wT) 



1/2 



exp{- 



2 sinh(wT) 



[0*2 



cosh((jjT) — 2x1X2]} 



(3.29) 



(3.30) 



is the euclidean propagator 119] for a particle of mass m moving, in a time interval T, under the harmonic potential 



2 muj x . The initial and final coordinates are x\ and x%, respectively. We obtain from (3.29) and (3.30) the asymptotic 
result 
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Z(X) 



da- 



(D A)i/2 



•exp( 



n 2 D Q R A A 2 
2L 2 ' ~~a~ 



1 

X 2 



(3.31) 



A simple way to understand the regularization of the divergent expression (2.4C) for Z{X) is that the additional 
terms in the path integral summation, associated to fluctuations, are complex quantities, which produce an increasing 
number of canceling factors as a — > oo. 



IV. SUBLEADING CORRECTIONS 



The asymptotic result ( |3.31 ) does not give us any dimensional parameter which could characterize in a more detailed 
way the circulation PDF, providing further motivation for a comparison with the experiment. We will investigate this 
problem here, through the analysis of subleading corrections for Z(X). 

Recalling what has been done, we observe that to derive expression ( p. 31 ) the path-integral for Z(X) has been 
written in a form which depends on an ordinary integral over a. The integrand is obtained from the saddle-point 
method, yielding a consistent result only in the time interval —1/a <t<0. In this way, fluctuations of the velocity 
field were com pletely neglected for t < —1/a (for t > they do not contribute to Z(X) due to causality). An improved 
form for (3.2E) may be found, thus, through the substitution 



G({x 2 \xi};T,uj,m) -> P(x 1 )G({x 2 \x 1 };T,uj,m,) , 



(4.1) 



where P(xi) is the probability density to have x = x\ at time t\ = —1/a. In other words, the effects of velocity 
fluctuations for t < —1/a are simply encoded in the PDFs for p and 4> a - It is important to note that these random 
variables are related to the circulation at different length scales. We may write, in fact, 



T R ee 2w 3 (t)TrR 2 = 2(p(t) + Mt))*R 2 



-.(«)- 



o ■ 



(4.2) 



Above, is the circulation evaluated for a circular loop of radius R in the xy plane, while refers in an analogous 
way to a loop of radius Rq in a plane perpendicular to the unit vector x a . These loops are centered at the origin of 
the coordinate system. From ( 3.30[ ) we see that as a — > oo only small fluctuations of cf> a and ui 3 become important. 
These fluctuations are associated to the core of the circulation PDF, which is modeled by a gaussian distribution, 



P(T r ) ~ exp(- 



A(r) 2 



(4.3) 



where "r" gives the length scale. This form of the circulation PDF for small r r is a phenomenolo gical ingredient in 
our analysis, well supported by numerical and real experiments [§,|o). Using (0-4^3) we rewrite ( 3.29 ) as 



Z(X) 



4>l 



(D A)V2 



f 3 ' 

da / dpdp Y[ d(j) a dcj) a exp(— - 



A 



47T 2 i? 4 



0§])G({&|fc}; 



2L 2 
L' 2 



a )(1 A(R) 2 ' 



\v 2 R\ 
A(i?o) 2 



^-)G({0 2 |0 2 };i,a. T7r 



)G({^|0 3 };-,2a,-^-) 
a 2L>o 



xG({p|p};i,2a,^(§) 4 ) 
a Uq R 



(4.4) 



In order to compute (4.4), a very convenient simplification of (3.30) follows from 



X EE X\e 2 — X2t 
— —i^E. 

x ee X\e 2 — x 2 e~ 



(4.5) 



which allows us to write 

G({x 2 \x 1 };T,uj,m) 



. UT J exp{-— — — - {x+ 2 + - (x~ 2 } 
2ir smh(wi ) J 2 smh(ti;J ) 2 2 



(4.6) 



It is also necessary to define W3 and 4> a in terms of p + , p ,0+ and <f> a . We have 
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^3 = P + 03 = 2s J h(2) [e 2 (P + + ^3 + ) - e \ P + 3 -)] , 
1 



. .Je0+ 2 -e-Vr, 2 ] ■ (4-7) 
2smh(l) 



Substituting (4.6) and (4.7) into (4.4), the gaussian integrals may be readily evaluated, giving 



> (4-8) 

where 

/3 ~ (16sinh(2))^A- 1 ~ 7.6A" 1 . (4.9) 



In the computation of (4.8) we have assumed that 

A(R )R 2 
A(R)R 2 

which is clearly verifi ed in practice g] 



> 1 , (4.10) 



We may interpret (4.8) as the asymptotic approximation to the lorentzian Z(X) ~ (A 2 + /3 2 ) -1 , which leads, on 
its turn, to the stretched exponential P(r) ~ exp(— /3|r|). The tail decaying parameter j3 is inversely proportional, 
therefore, to the width of the PDF's core, 2A. This agrees with Migdal's conj ectu re that P(r) is a function of the 
scaling variable r/A( 2fc_1 )/ 2fe , as discussed in the introduction. We would find ( |4.9| ) once again if we had considered 
other axisymmetric contours, as two concentric loops of radius R% and i?2, for instance. The PDF's dependence on 
the minimal area has to be completely contained in A, showing that universal features of the circulation PDF are 
related essentially to the form of its core. The manifestation of universality not only at the tails of PDFs seems to be 
in fact a property shared by other turbulent systems, as discussed recently in the problem of a passive scalar advected 
by a random velocity field in one dimension [|2l|| . 

A physical picture that may explain in more concrete terms the core-tail relationship for the circulation statistics, 
the result of the above computations, is in order. We may imagine that the large scale forces generate smooth 
configurations with small vorticity which are then fragmented in the cascade process up to the inertial range scales. 
These are the "soft" vortices that contribute to the core of the circulation PDF. With some probability, however, 
these vortices will be found in regions of the fluid characterized by high stretching. Their vorticity will be, thus, 
strongly enhanced, producing the intermittent configurations, described by the PDF tails. Since longitudinal velocity 
differences responsible for stretching do not fluctuate so quickly as the transverse ones related to circulation, the 
correlations of the soft vortices are transposed to a different range of vorticity. This is the meaning of /3 ~ A -1 , which 
implies that the same anomalous exponents determine the tails and the core of the circulation PDF. 

It is clear, from the results just obtained, that our task, within the reach of the saddle-point method, is at best to 
establish predictions suitable to experimental test, even if we lack a precise knowledge of A(i?), to which further and 
complementary investigations have to be directed. One might suppose that A(R) could be derived, at the onset of 
turbulence, from the viscous limit of the Navier-Stokes equations, in such way that the circulation PDF would keep 
the form of its core, while developing slowly decaying tails. In the viscous case, the circulation PDF is indeed gaussian, 
but A(R) ~ R 2 (see appendix C), which is in strong disagreement with observations. Thus, we do not expect smooth 
configurations of the velocity field to play any role in determining the core of the circulation PDF, even in situations 
close to critical Reynolds numbers. 



V. PARITY BREAKING EFFECTS 



Let us study now possible asymmetries between the left and right tails of the circulation PDF, caused by parity 
breaking external conditions. We will investigate here two simple models (which will be denoted henceforth by A 
and B, respectively): a fluid in rotation with constant angular velocity Q = luz and a fluid stirred by the force 
f a (x,t) = f a (x,t) + f a {x), where only f a (x,t) is random, being defined by (2.2). The static component f a (x) is the 
one responsible for parity breaking effects. In these models we will assume that the core of the circulation PDF is 
given by a shifted gaussian distribution, 
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P(r) ~ exp(- 



(r-r ) s 

A 2 



(5.1) 



with To <C A, and A being the same as in the situation where parity breaking co nditi ons are removed ( uj = f a (x) = 0). 
To simplify the notation, we took out the scale dependence of T, Tq and A in (5.1). 
Model A 

A turbulent rotating fluid, with angular velocity uj = ujz is described by a slightly different version of the Navier- 
Stokes equations (PJ), which takes into account the presence of non-inertial effects: 



dtv a + vpdpVa - 2uje 3a7 Vj - u 2 x^ = -d a P + vd 2 v a + f a 



(5.2) 



The c entrifugal force uj 2 xj^ may be absorbed by the pressure term. Following all the steps carried in sec. II, equation 
(2.36) becomes now 



b + 2ab - 2auj = -2ttL> A ( - ] e 2at 6(-t) 



(5.3) 



which is solved by 



b(t) =lu- 



irD \ ( R 



2a 



-) e- 2 °!' 



(5.4) 



while equation ( 2.23| ) still yields the same solution for v a (x,t), given by (2.32) (this is also true for model B; the 
distinction between the models is due only to different solutions for b{t)). Using ( |5.l| ) and (5.4), we obtain the 
corrected form of (4.4), which gives, after computations are done, 



1 



Z(A) ,x P! -;A,^)^(1 - exp(-2^|)^) 



(5.5) 



We find immedi ately from (5J3) the shift Y — > T + uj in the circulation PDF, as expected on physical grounds. Another 
consequence of (5.5) is that the tail decaying parameter (3 gets multiplied by a factor which is related to the shift Tq 
at the core of the circulation PDF. As To increases, the PDF tails become broader, apart from the overall shift by uj. 
Model B 

Expanding the static part of f a (x, t) in a power series around x = 0, we will have, up to first order, 

f a (x) = f a (0) + d [fi f a] xp + d {(3 f a} X{] , (5.6) 

where 



d{/3fa] = Tjidpfa + dafp) 



x=0 



x=0 



(5.7) 



The above expansion is physically associated to parity breaking mechanisms defined in the integral scales. As a 
conjecture, we expect that the induced modification on the instanton solutions will lead to a model independent 
description of parity breaking effects at the PDF tails. 

Let us consider here the case where <9[/3/ a ] = £30/3/0 > to get equations which are still invariant under rotations 
around the z axis. The strength of parity symmetry breaking is given by the external parameter / . The first and 
third terms in the RHS of ( p3.6| ) are absorbed by the pressure in the Navier-Stokes equations. Similarly to the analysis 
of model A, we write the equation for b(t), 



R 

b + 2ab= -27tD A ( — 



e zat 8{-t) + f , 



(5.8) 



which solution is 



w 2a 2a \L 



(5.9) 
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From this we obtain, instead of (2.39), 



s<°> = 



where the tt/2 rotation A — > iX was taken into account, and we have 

foL 2 



ttD K 2 



(5.10) 



(5.11) 



The result ( j5~T(i| ) may be quickly derived if we note that the only implication of (|i]) is the shift r — > T + irR 2 fo/a 
in the MSR action, lead ing t o S (0) -> 5 (0) + iXirr 2 f /a. 

Using now (bj) and (5.10) to correct (4.4), we get, through a direct computation, 



Z(X) 



(A + + (i 2 



exp(-2^|) 



/? 2 



A 2 [(A + i(3) 2 + /3 2 



(5.12) 



From the above expression for Z(X) we find that the right and left tails of the circulation PDF are described by 
P+(T) ~ exp(-/3+|r|) and P_(r) - exp(-/3_ |r|), respectively, with 



/?+ =/3+[exp(-2^)/3 2 + ^ 2 ]5 , 

/?- = -/3+[cxp(-2^)/3 2 + ^ 2 ]i . 
It is interesting to note that the product of the tail decaying parameters is approximately constant: 

0+ • /?_ = exp(-2^|)/? 2 ~ /3 2 . 



(5.13) 



(5.14) 



There is a compensation effect between the left and right tails, as the parity breaking parameter /q is varied. 



VI. CONCLUSION 



The problem of circulation statistics in fully developed turbulence was investigated through the Martin- Siggia-Rose 
formalism. An infinite set of axisymmetric instanton solutions follows from the saddle-point equations, which are 
labelled the component a zz of the strain field, a partially annealed variable. In physical terms, this means that the 
non-diagonal components of the strain tensor, related to circulation, are in fact the random variables which fluctuate 
against the quasi-static background defined by a zz . The asymptotic behavior of Z(X) =< exp(iAr) >, as well as 
its subleading correction, were found, leading to a stretched exponential description of the tails of the circulation 
PDF, a result in agreement with observational data. The core and the tails of the circulation PDF were seen to be 
intrinsically related. We estimate the tail decaying parameter (3 to be approximately equal to 7.6A -1 , with 2A being 
the width of the PDF's core. The numerical value in this estimate is related to the transition at time t ~ — 1/a 
between the saddle-point dominated regime and the free turbulent description of the fluid in the MSR formalism, 
which corresponds to have A = in ( |2.6| ). More generically, if the transition occurs at time t ~ —g/a, where g may 
be regarded as an adjustable phenomenological parameter, then we will have (3 ~ 4sinh(2g)2 A -1 . The relationship 
between /3 and A implies that universal features of the circulation statistics are determined essentially by the PDF's 
core, which, however, cannot be approached by means of the instanton technique. 

Parity breaking effects were also studied, as the ones which occur in rotating systems or in fluids stirred by parity 
breaking external forces. Well-defined predictions were derived, which we believe are within the reach of present 
numerical techniques, like the method of direct numerical simulations. 

On the theoretical side, the important problem to be addressed in future investigations is just the study of the core 
of the circulation PDF. It is likely that some explicit characterization of vorticity filaments will be necessary in order 
to study matters as anomalous exponents associated to intermittency and the minimal area conjecture. 
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APPENDIX A: TIME-DEPENDENT TRANSLATIONS 



The MSR action S(X), eq. (|jj), is invariant under the group T of time-dependent translations between coordinate 
systems, defined through 

x — > x = x — I dtu(t) , v a (x,t) — > v' a (x,t) = v a (x + / dtu(t),t) — u Q (i) , 
Jo Jo 

v a (x,t) — ► w^(x,i) = u a (x + / dtu(t),t) , Q(x,t) — ► Q'(x,t) = Q(x + / dtu(t),t) , 



P(x,t) -> P'(x,t)=P(£ + dtu{t),t) + u a (t)x a . (Al) 

J o 

We observe that T-symmetry holds in the MSR formalism whenever functionals of the velocity field are defined at a 
fixed instant of time, being also invariant under usual galilean tranformations = constant). 

Suppose we have a solution of the saddle-point equations with v a (x — 0,t) — ipa(t). A time-dependent translation 
may be applied to find another solution v' a (x,t) with v' a (x — 0,t) = 0, which yields the same saddle-point action. 
Our task is just to determine u(t) from 



v a { I dtu(t),t) = u a {t) 



A simple iterative procedure may be devised to find u(t). To start, we note that equation (A2) gives 

Ua(0) = V«(0) ■ 

Taking now the time derivative of (|A2|), we get 



u a (t) — up{t) dpv a (x + / dtu(t),t) — dt 1 v a ( / dtu(t),ti) 

J x=0 ^ 

At t = 0, we have, therefore, 

w a (0) - 1^(0) 9/9«a(^, 0)| £=0 - 9 t w Q (0,i)| t=0 = , 



= 



that is 



««(0) = ^(0) 9^(^0)1-0 + ^(0) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



We may proceed in the same way, considering expressions generated at each level of the iteration, to find time 
derivatives up to any order and use them to construct the Taylor expansion of u a (t) around t = 0. 



APPENDIX B: DESCRIPTION OF THE VORTEX SHEET 

We are taking fluctuations of Va\x,t) to have a discontinuity at the surface dfl, which encloses Q, a volume with 
typical size Rq. Note, in first place, that we may write 

v£ ] (x, t) = v< (x, t) [1 - F(x, f )] + v> {x, t)F(x, t) , (Bl) 

where v<(x,t) and v>(x,t) are given by (3.8) and ( 3.23| ), respectively, and 

Fix t) = i ^' ^ (B2) 

^ ' ' 1 0, otherwise. 

The idea now is to investigate the consequences of the incompressibility constraint, d a Va\x,t) = 0. This and (Bl) 
imply that 

d a v< (x, t) = d a v> (x, t) = , 

(v<(x,t)-v>(x,t))n a = . (B3) 



1G 



Above, n a = n ■ x a , where n is the unit norm al v ec tor p ointing outwards the surface (9f2. Writing n a = R a pxp /\x\ 
where R a p is a rotation matrix, we get, from (|3.8| ), ( 3.23 ) and (B2), 



x 1 R I [(a a (t) - a a (t)) + (b a p(t) - e a(T puj a (t))xp] = 



(B4) 



This gives a a (t) = a a (t) and R~^{b a p{t) — e acr puj a (t)) — M 1 p, where M = M{x) is an antisymmetric matrix. Since 
there is in any closed surface d£l at least one point where R a p = 6 a /3, we find that b a p(t) is also an antisymmetric 
matrix. Therefore, R a p is constant on dfl up to rotations around x, yielding n — x/\x\. To put it in another way, f2 
is a sphere of radius Rq. A convenient expression for b a p{t) is 



(B5) 



allowing us to define (3.24). 



APPENDIX C: CIRCULATION PDF IN THE VISCOUS LIMIT 



To study the viscous limit, we just neglect the convection term in the Navier-Stokes equations. As a re sult, we ge t 
an instructive example where the circulation PDF may be exactly found. The saddle-point equations (2.9) and ( 2.10| ) 
are now replaced by 



i(dtv a — vd 2 v a ) — J d 3 xD a fj(\x — x'\)vp{x' , t) 

i(d t v a + vd 2 v a ) = \e 3l3a ^6( ri _ - R)S(z)S(t) . 
rj_ 



(CI) 
(C2) 



The incompressibility constraints d a v a = d a v a = have also to be satisfied. Using (CI) and (C2), the saddle-point 
action in the MSR functional may be written as 



5(A) = -^jv-dx 



(C3) 



All we need to do, therefore, is to find v a (x±,z = 0, t = 0) = v a (x±,0). Applying (d t + vd 2 ) on equation (CI), we 
will have, integrating by parts and using (|C2|), 



[di-u 2 {d 2 Y)v a {x,t) = -F a (x,t) 



(C4) 



where 



F a (x,t) = -X I d?x'D a0 {\x-x' \)e^-^5{r' x - R)6(z')S(t) 
D \2irR 



l2 espaxpexpi-j^) 



In Fourier space, equation (C4) becomes 



(C5) 



(or +v z ie)v a {k,uS) = F a {k) 



(C6) 



We obtain, thus, 



v a (x, t) 
1 

47T^ 



(2tt) 5 



d kdu) „ a — — — exp(ifc ■ x + iuit) 



d A k 



F a (k) 
k 2 



exp(ifc • x ~ vk 2 \t\) 



Since we are interested to know v a (x±, 0), it follows, from (|C7j) , that 



(C7) 



v a (xx,0) = -. — / d 3 k °j - exp(zfcx • ^_l) 



k 2 



(C8) 
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Taking now (C5), we get 



F a {k)= / d°xe 3l3a x fj exp(-zfc • x - — ) 

, DqXtt^R 2 , L 2 k 2 , 
= -lespakp exp( — ) . 



(C9) 



Substituting this result in (C8), we will have 



nDnXR 2 



6u 



-£3f3aX{3 



Thus, from (C3) and (C1C), the saddle-point action is computed as 

X 2 D^ 2 R A 



5(A) = - 

Performing now the analytical mapping A — > iX, we find 



62/ 



X 2 D tt 2 R\ 

Z(X) oc exp( ) , 

w 

which leads to a gaussian statistics, described by the circulation PDF 

1 r 2 

where 



(CIO) 



(Cll) 



(C12) 



(C13) 



. 2Dn . i „n 
A = (^)iirR 2 



(C14) 
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FIG. 1. 

The three axisymmetric surfaces of revolution, I, II and III, which bound the support of small scale velocity 
fluctuations determined by S^ x > . As a — > oo, the surfaces asymptotically approach the cone given by z 2 = (x 2 + y 2 )/8. 
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